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Abstract 

A large class of polarized and unpolarized deep inelastic data is successfully described 
with Fermi-Dirac functions for the non-diffractive part of quark parton distributions. The 
NLO approach used here improves the agreement with experiment of the previous LO work. 
We get a broader distribution for the strange parton s{x) than for s{x). 



PACS numbers: 13.60.-r, 13.60.Hb, 14.20.Dli 



1 Introduction 



The possibility of having s{x) ^ s{x) has been advocated to explain the conflict between 
two different determinations of the strange quark sea in the nucleon. The CTEQ global analysis 
obtained the strange quark distribution on the assumption of s{x) = s{x). Actually, the 
measured quantity was the following combination of nucleon structure functions on an isoscalar 
target, A^: 

A (f^ + - 3 = f [six) + -six)], (1) 

where the sum of F2 neutrino structure function came from CCFR DIS data from neutrino 
and antineutrino beams Q and from the NM Collaboration Q. The result of this mean 
strange-antistrange quark distribution was found to be quite different from the strange quark 
distribution extracted from dimuon events of CCFR ||5| . Since the neutrino events dominate the 
CCFR data set, one should consider the latter result as s{x). Therefore one may consider the 
conflict as a first evidence for the difference between s{x) and six). Note that as a result of 
the reanalysis by CCFR including the effect of higher-order QCD corrections, the observed 
discrepancy is reduced, but it is still significant. An interesting contribution to this problem 
would be represented by an independent evidence of ^ six). New neutrino and antineutrino 
data on F2 and F3 have been measured and they are quantities (especially -F3) sensitive to 
six) 7^ six). Therefore a global analysis of nucleon structure functions, including F2 and -F3 
from neutrino and antineutrino beams rather than the dimuon events of CCFR measurement, 
may provide this kind of independent confirmation. 

On the theoretical side, with the heuristic argument that m\/3 > ttik/'^ Q a broader 
distribution for (expected to combine with the valence u and d quarks to give a A) than 
for s(x) (expected to combine with the same quarks to give a K) was advocated. This is just 
the property that could explain the previously mentioned data conflict. 

In another scenario, a different shape for six) and s(x) may be naturally understood in 
the framework of an approach developed in the last years ^, |^]. With the motivation of 
keeping into account the role played by Pauli principle in explaining several experimental facts, 
the parton distributions were described in terms of the sum of a gas component, given by a 
Fermi-Dirac function for quarks ii = u\ , d\ , u, d), 

A (1 - x)^ 

Pi,gas\X) — /„ -\ ~^ (2) 

exp [^) + 1 
and a liquid term, unpolarized and isoscalar, given by 

L(x) = 0.12x~i-i9 (l-x)9-6. (3) 
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For the gluons a Bose-Einstein function was taken {i = g\ g^), 

8 Ax" (1-x)^ 

Pi,gas[X) — , . ■ (4j 

3 exp - 1 

In Eq.s and @) x plays the role of a temperature, Xi of a thermodynamic potential, depending 
on the flavour and the spin of each parton, and the function in the numerator is a weight 
function for the density levels of the quarks in the Pz = oo frame of reference. In Ref. Q this 
statistical parametrization was applied for describing both polarized and unpolarized data and 
compared with a standard one, obtaining a very good agreement between the first, second, and 
third moment of the distributions. 

However, it is rather natural to assume different weight functions for q and q, since in 
the nucleon, which is not a C invariant object, the quarks (transforming as a 3 under SU (3)c) 
dominate over the antiquark (transforming as a 3). Indeed, inspired by the previously mentioned 
heuristic argument {m\/3 > m;^/2), we expect a broader weight function for the quarks than 
for the antiquarks. The strange quark is the best one to test this idea, since it is the lightest 
not valence parton. 

With respect to the analysis of Ref. Q, that was performed at Leading Order (LO) in the 
strong coupling, a^, here we include the Next to Leading Order (NLO) corrections. We study 
a large class of polarized and unpolarized data within the framework just described, with the 
important modification of taking different weight functions for quarks and antiquarks. 

The paper is organized as follows. In Section 2 the input parametrizations of the parton 
distributions are presented. Section 3 is devoted to the description of the experimental data 
used in our analysis and of the method of resolution of evolution equations, and Section 4 
reports our results and conclusions. 



2 Parton Parametrizations 

According to the hypothesis of a role of Pauli principle for quark parton distributions, we take, 
at Qq = 3GeV^, Fermi-Dirac functions for u and d quarks {i = u\ u^, d\ d^), 

Ql) = + I hix), (5) 

exp (^) + 1 2 

where the last term, 

/l(x) = Al x"Ml-^)^^ (6) 

represents an unpolarized and isoscalar part which takes into account the large diffractive con- 
tribution in the small x region (—2 < < — !)• 
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With the aim of ahowing different shapes for s{x) and s{x), we take a different form for the 
weight function of antiquarks. Since at small x the single parton contribution is overwhelmed 
by the diffractive one, it is very difficult to determine the difference in the exponents of x. 
So, we assume them equal for quarks and antiquarks, allowing, instead, different values for the 
exponents of (1 — x), which is sensitive to the large x behaviour where the diffractive contribution 
becomes negligible. In conclusion we have 

f,ix) = A.x'^ (l-xf^, (7) 
fgix) = Agx^ {l-xf\ (8) 

with a > —1. 

In the analysis of Ref. we found rather negative values for the potentials of the partons in 
the sea, which correspond to the Boltzmann limit for their distributions. We then parametrize 
s, u, and d in the following way: 

u{x,Ql) = ku fq{x) + fiix), (9) 
d{x,Ql) = h Mx) e-^ + hix), (10) 

s{x,Ql) = Kfg{x)e-^^+^fL{x). (11) 

Note that the different coefficient in front of /l(2;) in Eq. (^) takes into account the results of 
Ref. in the relative size of the strange sea with respect to the non-strange one. This result 
is also used to decrease the number of parameters, in writing 

s-(x,Qg) = ^^^l^i^±^^, (12) 

with, however, the requirement that the first moments of s and s be equal, since the nucleons 
have no strangeness. 

The defect in the Gottfried sum rule [0] implies a non trivial sea in the nucleons, with d > u, 
as confirmed also by the experiments NA51 at CERN |l2| and E866 at FNAL |l|]. Along this 
line, we should also consider the possibility of polarization in the quark sea, as advocated to 
account for the large defect in the Ellis and Jaffe sum rule |0] for the polarized structure function 
of the proton, g^, first shown in the EMC experiment ||l^ and confirmed by the following ones. 
We expect the distributions to be proportional to the gas component of the unpolarized ones, 

Au{x,Ql) = Aku[uix,Ql)- fLix)], (13) 
Ad{x,Ql) = Ak^[d{x,Ql)-h{x)], (14) 

As{x,Ql) + As{x,Ql) = Aks 



six, Ql) + s{x, Ql) - -^fUx) 



(15) 



with |AA;j| < 1 and for simplicity we have taken the same proportionality constant for s and s. 
The isovector combination Au — Ad is thus proportional to the gas component of -u + c? by a 
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factor, which is less than one in modulus for the positivity constraints on polarized distributions, 
and reaches that value only in the extreme cases of full and opposite polarization for u and d. 
For the isoscalar combination we have a part proportional to fq, Au + Ad and As, and another 
one. As, proportional to fq. In conclusion, we have three parameters, Aku, A/cj, and Akgg to 
describe the sea contribution to the two polarized structure functions of the nucleons, and 

For the gluon distribution we take the Bose-Einstein form {i = g\ g^), 

Mx,Ql) = /f^^ , (16) 

exp(2^)-l 

with fg given by 

fg(x) = ^Agx"^ (l-x)^^ (17) 

but with the constraint ag < a. Note that a more divergent term, similar to the second one in 
the r.h.s. of Eq. (^), seems to be excluded by data; so, we chose not to include it in Eq. (16). 



We impose that the total momentum of the partons be unity. Other constraints on the 
distributions concern the values of the following unpolarized sum rules: 

lAdler = 1.01 ±0.20, (18) 

Igls = 2.5 ±0.08, (19) 
iGottfried = 0.235 ±0.026. (20) 

As far as the Adler sum rule jl^ is concerned, we consider this constraint as an exact one, while 
we use the experimental errors for the other ones. Moreover, in calculating the theoretical value 



of the Gross-Llewellin-Smith sum rule [^] we include the QCD corrections up to O (a^) |18] 
iGLsiQ') = I'alsiQ' 



(21) 



where 



Mjq'^)= dx[u-u + d-d + s-s]{x,Q'^). (22) 
Jo 

Finally, we constrain the ratio {u/d){x = 0.18) to the value 0.51 ± 0.06 measured by the experi- 



ment NA51 112 



3 Description of Data and Parton Evolution 



We perform a NLO description, in the MS scheme, of the unpolarized data on F2 and F3 from 
CCFR 0, ¥2''^ from NMC ||l^, and of the polarized structure functions measured at SLAG in 
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the E142 E143 0, and E154 (2|] experiments, at HERA in the HERMES experiment 
and at CERN in the SMC experiment |24|. All the unpolarized data were subjected to the cuts 
> 2GeV'^ and VF^ > lOGeV'^. We add a 1.5% uncertainty to the statistical errors of CCFR 
data since no global systematic errors are given. 



For solving the DGLAP evolution equations |^5[ we use the Jacobi polynomial method |2q| , 
which some of us already used in an analysis of polarized structure functions |27]. Here, we 
briefly recall the procedure. The given structure function is expressed in terms of a truncated 
series of Jacobi polynomials, 

N 

F{x,Q')=x-{l-xr Y: a^f \Q')&f{x), (23) 



k=0 



where Q?^ satisfy the orthogonality condition 







x^il - xf efef^ dx = 5ki. (24) 



By using the definition of the coefficients, ci^^^\ 

a^:^\Q')= t F{x,Q')Qf{x)dx, (25) 
Jo 

with a little algebra it is possible to obtain for F{x, Q^) the following expression, 

N k 

F{x,Q^)= x'^il-xf J2 E (a,/3)F,+i(g2), (26) 



where Cj'^ (a, /?) are the coefficients of the expansion of the Jacobi polynomials, 0^^(x) , in power 
of x, and Fn{Q'^) are the Mellin moments of -F(x, Q^), 

FniQ^) = C x"^-^ F{x, Q2) dx. (27) 
Jo 

In this way, the dependence of F{x, Q^) is factorized in its moments, for which the solution 
of the evolution equations up to NLO is well known. 

In this analysis we reconstructed the unpolarized structure functions with = 12 and the 
polarized ones with A^ = 9. Moreover, we used different values of a and /3 for the different data 
sets. Table 1 reports the values of the parameters which give the best convergence of the Jacobi 
expansion. 

The unpolarized parton distributions at Qq are combined to give the following non-singlet 
(Qi), singlet (S) and gluon terms (G) (we suppress for brevity the x and Qq dependence of the 
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various terms), 



Qp — 

Qn — 

= 

Qs = 
s = 

G = 

while in the polarized case we have 

as 



2 {u + u) - {d + d) - (s + s), 
2 {d + d) - {u + u) - (s + s), 
u — u + d — d + s — 
s — s, 

u + u + d + d+ s + s, 
9, 



4 2 - 

SQp = a3 + ^ = - {Au + Au) - -{Ad + Ad + As + As), 

SQn = -as + ^ = ^ (Ad + Ad) - ^(Ati + An + As + As) 
5S = ao = Aii + Att + Aii + Ad + As + As, 
6G = Ag. 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 

(34) 

(35) 
(36) 
(37) 



The moments of the previous quantities at Qq are evolved to the of data (see [^] for the 
complete expressions of the evolved moments). The relation between the evolved moments of 
the combinations in Eq.s (p^)-(|37|) and the structure function moments, to be used in Eq. (p6|), 
comes from the expression of the unpolarized and polarized structure functions. Neglecting for 
the moment the charm contribution, we have for the first ones (for neutrino structure functions 
we give the expressions for isoscalar targets) 
Ff(x,Q2) 



X 

F2^^^«(x,g2) 



w . 



C| Q Qp + ^ I (x, Q2) + (c| ® G)(x, 
C| Q Q„ + ^ I (x, + (cf G)(x, g2), 



- {F^^{x,Q^)+Fl^{x,Q^)), 



(38) 
(39) 
(40) 



1 

X 



+ (2 a - 1) F!^ 



(v-v) 



1 + lK.I 
6 



(2 S + Qp + Q„)+ 



(S - Qp - Q„) + (2 a - 1) 



w . 



(Q3-Qs) + IK.I' Q. 



(x,Q') + 



(C|®G)(x,q2 



x,Q') 



1 + iKd 



(Q3-Q.) + |K.|' Q. 



(x,Q'). 



(41) 
(42) 



In the previous equations and Cf are the coefficient functions for quarks and gluons, which 
at LO are given by 

Cf(x)=5(l-x), C|(x) = 0, (43) 
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and at NLO can be found, for example, in 

'•I dz 



The convolution is defined as 



-/ 



9{z) 



1 dz 



(44) 



In the neutrino structure function the Cabibbo-Kobayashi-Maskawa |2£] elements appear, for 
which we use the following values: |KdP = 0.9508, = 1 - |KdP, iKdP = 0.0488, ^ 

0.9493. Moreover, a = 0.828 is the fraction of v with respect to v in the CCFR experiment l30|] 



and a term with a in the r.h.s. of Eq. (42) does not appear since the CCFR data have already 
been corrected for the s + s contribution P] . The polarized structure functions are 



1 



12 



(^C" 6Qp^n) (x, g2) + ® <5E) (x, 



5f(x) 



- [1--u^d] K(x)+<(x)] 



(45) 
(46) 



where (5C* are the polarized coefficient functions for quarks and gluons (see, for example, [31] for 
their expression) and ojd = 0.058 [^] is the D-wave component in the deuteron ground state. 
Note, however, that, when experimentally available, we fit the asymmetries 



2x (l + i?P'"(x,g")), 



(47) 



with R = Fl/{2xFi 



As stressed in |3^, a consistent treatment of heavy flavours can be carried out in the Fixed 
Flavour Scheme (FFS), where the heavy quarks are not considered as intrinsic partons, but 
produced by the interactions of the other partons (light quarks and gluons). To this aim, we fix 
the number of active flavours in the splitting functions to be 3 and add the charm contributions 
to neutral and charged current F2 and x F^ (we neglect the small h quark contribution and the 
c contribution to the polarized structure function (71), 



7(c) 

2 NC 
■ 2 CC 



(x,Q2 
(x,Q2 



— " X / C' 

TT 



F!t\x,Q') 



2. '(6. 1 _E n-f(f,QO^(.). 



(48) 
(49) 

(50) 



p=q ,a 



In the previous equations the hat indicates that we are calculating a quantity at = /i^ and 

C{x,Q^) = C{x,Q^,fi^), Hix,Q^) = H{x,Q^,fi^), (51) 

where /U is the factorization scale, equal to 4 and + in the neutral and charged current 
processes, respectively; the expressions of the C and H^s can be found, for example, in |34] and 
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5|; S,{x,Q ) = X a{Q^), where 



a(Q2) 



1 + i^, for NC, 



mi 



(52) 



l + -f, forCC. 



Furthermore, in Eq.s (|49|) and (|50| ) we have 



1 = \ 



(2 S + Qp + Q„) + -Qp- Qn) 

b 3 



(Q3-Qs) + |K.|'Q. 



(Q3-Qs) + |K.|'Q. 



for F^+^ 



for Fr", 



(53) 



for F 



3- 



Differently from the charm quark treatment, in the expression of the QCD couphng constant, 
we include the usual active flavours, nj, below each threshold. 



<"°(Q^) 

47r 

/3o = 11 - |n/. 



1 [3i Inln^ 

(3i = 102 - f n/, 



(54) 



where A^lo is fixed to the value 0.2263, so to have as(M|) = 0.118 [||]. Whenever IS necessary 
we use the following values for the heavy quark masses: 



rric = 1.5 GeV, rub = 4.5 GeV. 



(55) 



4 Results and Conclusions 



We consider the three following cases: 
Fit 1: Au, Ad, As + As / 
Fit 2: Au, Ad / 0, As + As = 
Fit 3: Au = Ad = 0, As + As / 

In Table 2 the values of the parameters for Fit 1 are reported. We do not show the results 
for Fits 2 and 3 because they look very similar to Fit 1 in the parameter values and Xred (1-76 
and 1.80, respectively), with a very small difference in the gluon contribution, which results in 
a negligible positive polarization (Ag = 0.058) for Fit 2 and a slightly larger one {Ag = 0.113) 
for Fit 3. Table 3 shows the comparison between the polarization of antiquarks in Fit 1, 2, 
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and 3. Interestingly enough, u in Fit 2 and s + s in Fit 3 come out fully negative polarized 
to confirm that, as supposed in the interpretation of the EMC result in the gauge- invariant 
factorization schemes, the sea is negatively polarized. It is however difficult to disentangle 
the sea contribution to the polarization from the valence one, expecially since sea partons are 
mostly present at small x, where the large diffractive contribution makes more unprecise the 
determination of the polarized structure functions. Indeed, one has to find the polarized cross- 
section, which is expected to vanish in the limit x — > 0, from the difference of two cross-sections, 
which go to infinity in the same limit. The values of = An + Au + Ad + Ad — 2{As + As), 
which are related from SU{3) symmetry to the combination 3F — D = 0.579 =b 0.025 are 
very different for Fit 2 (0.373) and 3 (0.884), while Fit 1 is similar to Fit 3 (0.956). This shows 
that at the moment we cannot provide a test for this SU{3) prediction. 



In Table 3 the values of the polarized sum rules, Bjorken [38|, Ellis- Jaffe [14| for the proton 
and for the neutron, are also reported in correspondence of the three considered cases. The values 
for the Bjorken sum rule are consistent with the 0{as^) theoretical prediction at = 3GeV^, 
0.174 ±0.002. 

The quark distributions at Qq are plotted in Fig. 1. In Fig. 2 we compare Qp, Qn, S, and 
g, evolved to = IGeV^, with the same quantities obtained from the global fit of Ref. ||39t| . 
Note, however, that in Ref. [^9| a different prescription for the charm quark treatment, 
instead of the FFS, was usedQ and the value rric = 1.35 GeV. Fig.s 3-13 show the comparison 
with experimental data of the structure functions calculated for the values of parameters of Fit 
1 (solid lines). We get a very good description of polarized data, which, however, have larger 
errors than the unpolarized ones. As far as the latter are concerned, a similar good quality is 
exhibited by the neutrino structure function x F^, while the description of F2 is not so accurate 
at low X. To some extent this can be observed for the high-x description of F2 and F2 too. Note, 
however, that if the CCFR data are not included in the analysis, the high-x behaviour of F2 
and F2 drastically improves {Xred ~ ^-^^ for the values of parameters reported in Table 2 as Fit 
la), as is evident from the comparison with data of the dashed lines in the figures. On one side, 
this suggests that, even with our hypothesis on strange quarks, the discrepancy between NMC 
neutral current data and CCFR charge current ones remains unresolved, and one has to consider 
other effects like charge asymmetry ||l], |3^. On the other side, the high-x difference between 
solid and dashed lines in Fig. 3 and 4 is also justified by the fact that if the experimental ratio 
between i^J ^-^id x F3 at high-x is different from 4/9 (the value implied by the dominance of 

in that region) one cannot describe in a fully satisfactory way both the structure functions 
unless allowing more freedom in the quark parametrizations. 

With respect to the results found in our previous LO analysis P] we can confirm the pattern 



'Actually, in Ref. a variant of the ACOT scheme |4l| was implemented. See Ref. for an extensive 
review on this issue. 



10 



of the ratios "p^^ between the second and the first moments of the gas component of the 
distributions of n^^, n^, S ^ and d\ quarks, and we get similar values for the temperature^ x. 
Moreover, in agreement with the conclusions by Brodsky and Ma we find a broader s 
distribution than s, as shown by the higher value of the second moment of s with respect to s. 
More precisely, we have for Fit 1 (similar values result for Fit 2 and 3) 

+ I 0^0421 for 

4.2 \ 0.0395 for », ' ' 

is the second moment of the liquid component) with a relative difference of ~ 7%, compared 
to the value of 10% obtained in [|l|. 

Finally, we can conclude that the fact that with the present NLO analysis we find Xred ^ 
(and < 1 for Fit la), which is smaller than the value found in our previous approach is a 
good point in favour of statistical distributions. 
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a 


/3 


NMC 


-0.078 


0.077 


FcCFR 


-0.99 


0.500 




-0.096 


1.210 


pec 
c 


-0.096 


1.210 


Fi 


-0.503 


0.006 


91 


1.084 


3.550 



Tabic 1: Values of a and /? in the Jacobi description of the structure functions relative to the 
NMC and CCFR experiments, to the charm component in neutral and charge current processes, 
and to Fi and gi, respectively. 





Fit 1 


Fit la 


^red 


1.76 


0.97 




3.55 (1 _ a:)2-76 


3.34 (1 - a;)2-62 


fq{x) 


7.33 a;-0-275 (i _ a;)4-84 


8.28 a;-0-254 (I - x)^^-^ 




87.4 a;-0-276 (i _ a;)7-36 


88.0 a;-0-255 (i _ a;)8-57 




0.124 (l_^)14.7 


0.120 (l-a;)i2-6 


X 


0.250 


0.257 


AG 


-0.009 


-0.193 








P?'/Pi'' 


Xi 








0.750 


1.58 


0.149 


0.867 


1.50 


0.161 




0.150 


0.870 


0.118 


0.184 


0.836 


0.126 




-0.033 


0.581 


0.109 


0.047 


0.639 


0.119 




-0.136 


0.440 


0.10-5 


-0.100 


0.447 


0.113 


a' 


-0.427 


3.88 


0.055 


-0.409 


3.93 


0.052 


gl 


-0.426 


3.89 


0.055 


-0.399 


4.12 


0.052 




ki 






ki 






U 


0.099 


0.181 


0.073 


0.118 


0.155 


0.045 


d 


0.181 


0.330 


0.073 


0.232 


0.303 


0.045 


3 


0.114 


0.122 


0.095 


0.112 


0.109 


0.101 


S 




0.122 


0.073 




0.109 


0.045 



Table 2: Values of the parameters of the input distributions for Fit 1 and Fit la (see text). We 
denote the first and second moments of the gas component of the distributions with pf'^ and 
p\^\ respectively. 





Fit 1 


Fit 2 


Fit 3 




Aki 




(2) 
Pi 


Ak, 




(2) 
Pi 


Aki 




(2) 
Pi 


Au 


-1.00 


-0.181 


-0.013 


-1.00 


-0.185 


-0.013 








Ad 


0.24 


0.080 


0.006 


0.50 


0.017 


0.001 








As + As 


-1.00 


-0.243 


-0.020 








-1.00 


-0.208 


-0.018 


as 


0.956 


0.373 


0.884 


Bj 


0.174 


0.175 


0.176 


EJp 


0.133 


0.134 


0.136 


EJn 


-0.041 


-0.041 


-0.040 



Table 3: Values of the polarization parameters and sum rules for the three cases considered in 
the analysis. 




Figure 2: Comparison of Qp (solid lines), Q„ (dashed lines), S (dotted lines), and g (dash-dotted 
lines), evolved to = iGeV'^ (BPR), with the same quantities obtained from the fit of Ref. 
H (MRST). 




Figure 3: Comparison of the prediction of the fit with the experimental data on F2 from NMC 
[19|. For display purposes -Fg been multiplied by the numbers in brackets. The solid and 
dashed lines correspond to the value of parameters for Fit 1 and la, respectively (see text). This 
notation is hereafter adopted. 




ure 4: Comparison of the prediction of the fit with the experimental data on F2 from NMC 
. For display purposes F2 has been multiplied by the numbers in brackets. 




ure 5: Comparison of the prediction of the fit with the experimental data on F2 from 
For display purposes the numbers in brackets have been added to F^'. 
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Figure 6: Comparison of the prediction of the fit with the experimental data on x from CCFR 
[^. For display purposes the numbers in brackets have been added to x F^. 




Figure 8: Comparison of the prediction of the fit with the experimental data on Af from E143 
The Unes are evaluated at the of the experimental points. 




Figure 9: Comparison of the prediction of the fit with the experimental data on 5" from E154 
[p2| . The Unes are evaluated at the of the experimental points. 




Figure 10: Comparison of the prediction of the fit with the experimental data on from E154 
[p2| . The lines are evaluated at the of the experimental points. 




Figure 11: Comparison of the prediction of the fit with the experimental data on from 
HERMES The lines are evaluated at the of the experimental points. 
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Figure 12: Comparison of the prediction of the fit with the experimental data on from SMC 
p4[ . For display purposes data at adjacent x values have been grouped together and the numbers 
in brackets have been added to vl?. 
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Fi gure 13: Comparison of the prediction of the fit with the experimental data on A'^ from SMC 
p4[ . For display purposes data at adjacent x values have been grouped together and the numbers 
in brackets have been added to A'f. 



